LU factorizations, q = limits, and p-adic interpretations of 
some g-hypergeometric orthogonal polynomials 

Tom H. Koornwinder and Uri Onn * 
Dedicated to Dick Askey on the occasion of his seventieth birthday 



Abstract 

For little g-Jacobi polynomials and g-Hahn polynomials we give particular g-hypergeometric 
series representations in which the termwise q = limit can be taken. When rewritten in 
matrix form, these series representations can be viewed as LU factorizations. We develop a 
general theory of LU factorizations related to complete systems of orthogonal polynomials 
with discrete orthogonality relations which admit a dual system of orthogonal polynomials. 
For the q — orthogonal limit functions we discuss interpretations on p-adic spaces. In the 
little O-Jacobi case we also discuss product formulas. 

1 Introduction 

This paper is concerned with limits for q J. of some g-hypergeometric orthogonal polynomials, 
in particular little g-Jacobi polynomials and g-Hahn polynomials. Limits of g-hypergeometric 
polynomials as q \ 1 are well-known, see [131 Chapter 5]. Many (<7-)hypergeometric orthogonal 
polynomials have interpretations in connection with (quantum) group representations, for in- 
stance as spherical or intertwining functions, matrix elements, Clebsch-Gordan coefficients and 
Racah coefficients, see for instance |29| and |12| . Often, the q j 1 limit of the polynomials 
corresponds to the q f 1 limit from the quantum group to the classical group. 

Limits of g-hypergeometric polynomials for q j have been considered for g-ultraspherical 
polynomials (see |TJ §5]) and for more general Askey- Wilson polynomials (see [U pp. 26-28] 
and references given there). The limit functions have interpretations as spherical functions on 
homogeneous trees (see references in ^ p. 28]) and on infinite distance-transitive graphs (see 
|3Uj ) . Note that homogeneous trees are locally compact but noncompact homogeneous spaces 
of the group GL(2, Q p ) (Q p the field of p-adic numbers). No geometric explanation of this q J, 
limit is known, see also the discussion in |14j . Macdonald considered the q j limit of Macdonald 
polynomials, both for root system A n (yielding Hall-Littlewood polynomials, see ^2 Ch. Ill]) 
and for general root systems (see |19| §10]). The limit functions have interpretations as spherical 
functions on a p-adic Lie group (see ^H)j m particular in the Ad-i case on GL(d,Q p ) (see [El 
Ch. V]). 
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Haran considered limits for q J, of g-Hahn polynomials, little g-Jacobi polynomials 
and little g-Laguerre polynomials. He gave interpretations of these limit functions as spherical 
or intertwining functions on compact p-adic groups; in particular in the little O-Jacobi case on 
the group GL(d,Z p ) (Z p the ring of p-adic integers). In Haran's ideas were extended to 
the higher rank case and the importance of the cellular basis, defined in [H], was emphasized. 
All these interpretations are valid with respect to the larger family of p-adic fields, or more 
generally non-Archimedean fields. Altogether, ^-special functions clearly play an important role 
as interpolants between representations of groups over all local fields (R and C for q = 1, and 
non- Archimedean local fields for q = 0). 

The present paper gives new proofs of Haran's 11 q [ limit results for little g- Jacob i 
polynomials and g-Hahn polynomials (sections [2] and 0J) by starting with a series representation 
for these polynomials which can be viewed as writing the square (possibly infinite) matrix 
corresponding to the orthogonal polynomials as a product of a lower triangular matrix and an 
upper triangular matrix (an LU factorization). The matrix elements of the upper triangular 
matrix (simply (q x ',q~ 1 )k) are the q- analogues of the cellular basis in the rank one case of [HJ. 
This observation was decisive for the second author in order to find the higher rank analogue 
of this q I limit, see |21| . For the limiting little 0-Jacobi functions we give a product formula. 
In section El we discuss the LU factorization more generally for orthogonal polynomials with 
discrete orthogonality relations which form a complete orthogonal system and for which the 
dual orthogonal system also consists of orthogonal polynomials. We consider an upper-lower 
factorization as well, When these systems of orthogonal polynomials are moreover finite (so- 
called Leonard pairs) then our theory is related to Terwilliger |27j . In section we apply the 
general theory of section |3] to the little g-Jacobi and the g-Hahn case. Finally, in section El we 
give p-adic group interpretations of the q = results obtained in section |2J 

In an earlier version of this paper we also discussed big g-Jacobi polynomials. However, 
this family does not fit nicely into the general theory of this paper, and interpretations on p-adic 
groups are yet missing. 

Acknowledgement We thank Erik Koelink, Michael Schlosser, Paul Terwilliger, Michael Voit 
and the referee for helpful comments. 

2 Little g-Jacobi polynomials 

In this section we describe the main themes of this paper through the little g-Jacobi polynomials. 
These polynomials and their p-adic limit were the pioneering example which motivated this 
paper. The cellular structure of the p-adic Hecke algebra of Grassmannians of lines (see ^SJ is 
essentially an LU factorization of the matrix corresponding to the spherical functions w.r.t. a 
geometric basis. As this basis is the limit of little g-Jacobi polynomials, it was natural to seek 
for such factorization for the matrix corresponding to the latter polynomials. 

Throughout we assume < q < 1. See standard formulas for little g-Jacobi polynomials in 
Q3J §3.12]. 
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2.1 Limit for q 1 

Little q-Jacobi polynomials are given by 

fq~ n ,abq n+1 \ 
p n (x;a,b;q) := 2 <M ' ;g,gx (n G Z> ). (2.1) 

V ag J 

For 6 = they are known as little q-Laguerre polynomials or Wall polynomials p n (x;a;q) := 
p n (x; a,0; q), see []2J §3.20]. It will turn out that we have to rescale the parameters a and b in 
order to be able to take the limit of these polynomials for q j 0. We define: 

pf«{x) -p^q-^q-^q) = 2 x ; g, g* +1 ) (2.2) 

(ra G Z> , x G Z> U {oo}, < a < 1, 6 < 1). 
By |131 (3.12.2)] the functions (|2.2j) satisfy the orthogonality relation 

OO g 

= (m,n G Z> ), (2.3) 

where 



(a6;g)oo (g;g)a, 
w a,fe;g . 1 - afrg 2 "' 1 (a, a6; g) n 

1 - a&g"" 1 a" (g, 6; q) n ' 1 ' J 

Note that the weights Wx ,b ' q and the dual weights w^'' are positive under the constraints for 
a and b given in ()2.2j) . Since the little g-Jacobi are orthogonal polynomials with respect to 
an orthogonality measure of bounded support, they form a complete orthogonal system in the 
I? space with respect to this measure, so the functions (|2.2|) also form a complete orthogonal 
system in £ 2 (Z> ; w°" h 'i). 

As was observed in gl (4.1), (4.2)], [TUl Remark 3.1] and 5, (5.1), (5.3)], little g-Jacobi 
polynomials can alternatively be expressed as a terminating 3^1 by application of the trans- 
formation formula (III. 8) or Exercise 1.15 (ii)] to the terminating 2</>i in ()2.2j) . We obtain: 

*•»(«) = „><-'> (-«)» ^ 3^ £l) . (2 . 6) 

[a;q)n \ b a ) 

As was also observed in the papers just quoted, this is related to the fact that little g-Jacobi 
polynomials are the duals of g _1 -Al-Salam Chihara polynomials 

Qn {\{aq~ x + a- l q x );a, b | g~ x ) 

= (-l)Vg-H"- 1 ) {{ab)- l ;q) n (^"'^^f ; g, g^aft" 1 ) . (2.7) 
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From ()2.6j) and (|2.7|) we get indeed the duality 



-a) q 



n -in(n-l) v Q ' gjw a,b;q 



(M)n " V ' (b;q) x 
xQ x ^(qabyHq' n +abq n - 1 y,(qab)- 1 2,(qab- 1 ) 1 2-,q- 1 ^ <V^>„) (2.8) 

Formula ()2.6|) can be rewritten as: 

min(n,a;) 



i(n-fc)(n-fc-l) £_ a ^n-ft (^g i g)n-fc (<%" j_ g)fc , 



k=0 



(a;q) n (q;q) k 



We obtain as an immediate corollary of (|2.9|) : 

Theorem 2.1. T/ie limit functions (little O-Jacobi functions) 



exist. They are equal to 

a,b;0 / \ 

Po \ x ) 



a,b;0 / \ 
Pi [X) 



p a /'°(x) :=limp^(x) (x£Z>„) 



1. 



o(l-6) 



ifx = 0, 



I- a 

1 ifx>0, 



(q n ,q x ;q x ) k [x E Z> ). 

(2.9) 

(2.10) 

(2.11) 
(2.12) 



o 

a 



i/ < x < n — 1, 
if x = n — 1, 



1 -a 

1 if x > n — 1 



(n > 2). 



(2.13) 



Remark 2.2. Theorem 12.11 was first stated by Haran [111 (7.3.37)], where the limit functions 
(|2.11|) - ()2.13|) are given in [111 (4.4.9)]. The limit result there follows from the expression for 
little g-Jacobi polynomials on [111 p. 59, second formula from below], which reads in our notation 
as: 



< A9 (z) 



(b;q)r 



- 1 (a-v~ n ;g)n 392 V b,0 

min(n ' x ) / \n—k (v k \ 

{a;q) n -k {q;q)k 



q n ,a ^-q 1 n ,q x 



(2.14) 



fc=0 



Formula (|2.14|) follows from (|2.2fl by the transformation formula 01 (1.5.6)]. Theorem 12.11 can 
be obtained from (|2. 15f) by letting q J, 0. 
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Haran's [111 Ch. 7] notation is connected with ours by: 

4>q,n (0) 



(2.16) 



Haran ^2 pp. 61-64] also considers the little g-Laguerre case (a — > oo in (|2.16|) : 6 = 
everywhere in our notation) and its q = limit. 

Remark 2.3. From (|2.9f) we also obtain the following asymptotics of Pn' b,q (x) as q [ 0: 



1 - a 



(l<x<n-l) 



(2.17) 



Alternatively, ()2.17j) can be derived from the (/-difference equation [131 (3.12.5)] by induction 
with respect to x, starting at x = 1. Theorem 12. II can also be proved by use of (|2.17|) . 



From (|2.4j) and (|2,5|) we get limits 



<' 6;0 := lim w a /' q 

qlO 



LO% b ''° := lim u>% b ' q 
?|0 



1 - a 
1 - ab 

(i-q)(i-ft) 

1 -ab 
1 

1 - a 



a{l-b) 
(1 - q)(l - ab) 
o n (l - b) 



if x = 0, 

a x if x > 0, 

if n = 0, 
if n = 1, 



(2.18) 



(2.19) 



if n > 1. 



Note that ?4' fe;0 > (x G Z> ) and cj£' fe;0 > (n € Z> ) if < a < 1 and 6 < 1. The 
orthogonality relation ()2.3|l remains valid for g = 0, as can be verified by use of (j2,ll|) - (|2.13j) 
and l|2.18|) — (|2.19|) . Formally, we can obtain the case q = of ()2.3|) by taking termwise limits. 



2.2 LU factorization 

Formula (|2.9|) has the big advantage over ()2.15j) that it can be rewritten in matrix form as a 
product of a lower and an upper triangular matrix: 



P q = L»U», i.e., Pl x = J2 L lkU q k , 



k=0 



min(n,x) 
k=0 



where 



I'?,, ■ iC 1 " 1 ^). L 



n,k ' 
1 



\(n-k){n-k-\) 



UL := {q X ;q- L ) k . 



J >0) 



(2.20) 



,n-k (bq k ;q) n -k(abq n 1 ;q) k u 



(a;q) n (q;q) k 



(q r 



(2.21) 
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As a limit case of (|2.2U|) . (12.211) formulas (|2.11|) — (|2.13|) can similarly be rewritten in matrix 
form Q2.2Uj) with q = 0, where 



n,x 
u k,x 





n,n— 1 



r0 

n,n 



V 



a,b;0 

1 





if < 
if > 

if k ^ n or n 
o(l - 6) 



1 - a 
a 



1 



1-ab 
1 - a 
1 

1 - a 



1. 



1, 



if n 

if n > 1, 

if n = 0, 
if n = 1, 

if 77, > 1. 



(2.22) 



Again, L° is a lower triangular and £7 an upper triangular matrix. The functions c q k {q x ) := 
U k,x/ U k,k = (<?\<r l )k/{<l\q)k and c° k (x) := U% x (see ((QP and C22» can be considered as 
forming a cellular basis in the terminology of |3 §3.3] and The functions q x — » const. £7? 
can also be considered as the one- variable cases of Okounkov's [SOJ shifted Macdonald polynomi- 
als. 



Remark 2.4. For < q < 1 we can consider P q as the matrix of a unitary operator from the 
Hilbert space ^ 2 (Z> ; (w q )~ 1 ) onto the Hilbert space £ 2 (Z> ; uj q ). Let v fc := a~ k . Probably, U q 
is the matrix of a bounded linear operator from £ 2 (Z>o; (ii; 9 ) -1 ) to £ 2 (Z>o;t>)> and L 9 is the 
matrix of a bounded linear operator from £ 2 (Z>o;t>) to ^ 2 (Z>o;w <? ). Probably, the operators 
corresponding to U q and L q have bounded inverses with matrices given by the explicit matrix 
inverses oiU q and L q . 

2.3 A product formula for little O-Jacobi functions 
Theorem 2.5. The functions pn b '° satisfy the product formula 

oo 

P a n b '°(x)p a n b ->°(y) = 53 c a x i° zP a /'°(z), (2.23) 

2 = 
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where c^ b y° z is given by 



-x,y,z 



Oz,min(x,y) 

l-2a + ab 
1-a 

(1 - b)a z 


1 - 2a 



I- a 

.z—x 



a 



ifx^y, 

if x = y = z = 0, 

if = x = y < z, 
if x = y> z, 

if x = y = z > 0, 

if0<x = y<z. 



In particular, 



P a n b ] °(x)p a n ' b '°(y) =K ,b; °(min(x,y)) ifx + y. 
Under the constraints < a < 1, 6 < 1 we have Cx ,b y° z > for all x,y,z € Z>o iff 

< a < \ and 2 - aT x < b < 1. 

Proof Straightforward verification by (|2,ll|) - (|2.13j) . 

Corollary 2.6. The functions pn' b, ° satisfy the product formula 

1 — ab 



(y) 



1 - a 



E^<' 6,0 K' 6;0 (-), 



2 = 



where 



/~ia,b,0 



min(a;,i/,z)+l 
1 — ° ",,6,0 



E K' b ' (^)Pn A0 (y)K A0 (^X 



a,fe,0 



n=0 



1 -a& 

is symmetric in x, y, z, and for x < y < z explicitly given by 

if0<x<y<z, 

1 if0 = x = y<z, 

(l-b)~ l a~ x if0<x = y<z, 
th 

if0 = x = y 



/^-Ya,6,0 

^x,y,z 



l-2a + ab 



1 - a 
1 - 2a 



I (l-a)(l-6) 



a x if0<x = y 



z, 



z. 



(2.24) 



(2.25) 

(2.26) 
□ 

(2.27) 
(2.28) 



(2.29) 



Remark 2.7. Dunkl and Ramirez [Sj obtained the little O-Laguerre functions Pn b, ° for (a, 6) = 
(p^ 1 , 0) as spherical functions on the ring of p-adic integers, and they derived the above product 
formula for those special parameter values from that interpretation as spherical functions. 
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3 LU factorizations: the general case 
3.1 Lower times upper 

Let us put the results of ^12. 21 in a more general framework. Let Af := {0, 1, . . . , N} or Z>o and 
let Y := {y x }x£Af be a countable subset of R. Let {p n }n&N be a complete system of orthogonal 
polynomials on Y with respect to positive weights w x on the points y x : 

^2 Pn(Vx) PmiVx) Wx = K)" 1 6 n ,m (n,m G Af), (3.1) 

xeAf 

where u n > for all n G Af and where, by completeness, we also have the dual orthogonality 
relation 

y~] Pn{Vx) PniVx') u n = (Wx)' 1 5 X>X > (x,x G Af). (3.2) 
neJ\f 

Note that completeness will certainly hold if Y is finite or bounded in R. 

We will now introduce some square matrices with row and column indices running over Af. 
Let P be the matrix with entries P UtX : = Pniyx) ( n , x G Af). Let also W be the diagonal matrix 
with diagonal entries w x (x G Af) and let f2 be the diagonal matrix with diagonal entries uj n 
(n G A^). Then Ij^.lj) and (|3.2[) can be written in matrix form as, respectively, 

PWP 1 = ST 1 , ptfLP = W~ x . (3.3) 

Define polynomials c k of degree k by 

k-l 

c k (y) ■= H(yj - y) (k e Af, y g R), (3.4) 
and let C be the matrix with entries given by 

k-l 

Ck, x ■= c k (y x ) = Y[{yj ~ Vx) (k,xeAf). (3.5) 

3=0 

Then C k , x = if k > x, so C is an upper triangular matrix. Then, for certain unique coefficients 
B ntk with B n , n / Owe have: 

n 

= Y. B ^ kCk ^ (rieAf,yeR), (3.6) 

fc=0 

min(n,x) 

Pn,x = ^ B nik C k:X (n,xeAf), (3.7) 

fc=0 

P = BC, (3.8) 
where -B is the lower triangular matrix corresponding to the coefficients B n ^ k (n > k). 
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Both B and C have two-sided inverses because they are triangular matrices with nonzero 
diagonal entries. Furthermore, if P = B'C is another factorization of P with B' lower triangular 
and C upper triangular, then B' = BD, C = D~ 1 C for some invertible diagonal matrix D. 

Things become even nicer if we know that there exist orthogonal polynomials dual to {p n }, 
i.e., if there exist polynomials r x of degree x (x £ M) and a subset {z n } ne j^ of R such that 

Pn(Vx) = r x {z n ) (n, x G j\f). (3.9) 

Then the polynomials r x {x € M) will form a complete system of orthogonal polynomials with 
respect to the weights tu n on the points z n (n G Af). Thus we can apply the previous result to 
the r x . Put 

k-l 

B n ,k ■= Y[(zi - z n ) (n,k£j\f). (3.10) 

2=0 

Then B is a lower triangular matrix and for some upper triangular matrix C 1 we have 

min(n,a;) 

Pn{yx) = r x {z n ) = B n;k C' hx (n,xeN). (3.11) 

fc=o 

Then P = BC'. Hence 

P = BDC (3.12) 
for some diagonal matrix D with nonzero diagonal entries 5k- Hence 

min(n,x) k~l k-l 

Vn{Vx) = ^2 S k~[]_(zi- z n ) Y\_{Vj -Vx)- (3.13) 

fc=0 i=0 j=0 

3.2 Inverting the matrices B and C and computing 5k 

The following theorem is the special case a,j = 1, bj = of the Theorem in [161 p. 48], and it 
is also the case f(x) := (xi — x) . . . (x m _i — x) (m < n) of [221 P-54, Exercise 97] (we thank 
Michael Schlosser for this reference), but we will give here an independent proof. 

Theorem 3.1. For distinct complex numbers y n (n € N) let C = (C m)n ) mtn ^j\f be an upper 
triangular matrix given by H3.5jl . Then 

n 

(c- 1 )^ = n - y^y 1 (° < k < n )- ( 3 - 14 ) 

3=0 

Proof Let m < n. Put V m ^ n := Yl m<i< j <n (yj — yi)- We have to show that 

n m—l n 

n ~ yfc ) n ~ y^" 1 = 5m > n - ( 3 - is ) 

k=m j=0 j=0 
3^k 
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This is clearly true for m = n. For m < n the left-hand side of (|3.15|) can be rewritten as 



E tlivj-vk)- 1 

k=m j=m 

n 



k—m 



k=m 



(^r'Ec-i; 



k—n 



k=m 



m<i<j<n 

1 

Vm 



.n—m—l 
am 



1 1 



n-m-l n-m-l 
iJk-l Uk+l 



1 



n—m—l 
an 



.n—m—l 



„. n—m—l 
an 



0. □ 



It follows by taking transpose in (|3.14|) that the inverse of the lower triangular matrix B 
given by (|3.10j) is the lower triangular matrix B~ l with 



(B- 1 )^ = J] ( z i - (0 < ^ < rn). 



(3.16) 



i=0 
iy^n 



Put y v := yN if N = {0,1, ... , N} and put y v := lim^oo y x if J\f = Z>q and if the (finite) 
limit y v exists and is not equal to any y x (x 6 J\f). We will now derive an explicit expression for 
the coefficients 5^ in (|3.13|) involving p n (yu) for all n G Af. This expression will have practical 
usage if a simple explicit expression for p n {y y ) is known, as is the case in most examples. 

It follows from ()3.13|) that 

n k—1 k— 1 n k—1 

Pn{Vv) = E^ 11^ ~~ Zn ) II ^ ~~ Vv ^> = ^Z B n,kh \\iVj ~ Vu)- 

k=0 i=0 j=0 k=0 j=0 



Hence, by matrix inversion, 



m—l 



E ( B l )m,nP n ^) = 6 ™ II " ^ 
n=0 j=0 

Thus, by CTTBft wg obtain th.6 following formula, for 6 m '. 



llf=0 Vu) n= o ,=0 



(3.17) 
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3.3 Upper times lower 

We obtain from (|3,3|) that P = Q~ 1 (P t )~ 1 W~ 1 and from (|3,12|) (only formally in the infinite 
dimensional case) that = (B*) -1 !) -1 ^) -1 . Hence, 

P = n- 1 (B t )- 1 D- 1 (C- 1 ) t W~ 1 . (3.18) 

From Theorem 13 . 1 1 (|3 ■ 1 ( j |> and ()3.5j) we see that 

k 

((^) _1 )n, fc = fl^-^r 1 (0<n<k), (3.19) 

i=0 

(C 1 " 1 )^ = n^-^)" 1 (0<x<fc). (3.20) 

3=0 

When we substitute everything in ()3.18[) then we obtain 

k k 

p n {y x )=u- 1 w- 1 ^n^-^IIfe-^ 1 ' ( 3 - 21 ) 

fc>max(n,i) i=0 j=0 

Remark 3.2. The derivation of (|3.18|) is purely formal if J\f = Z>o, since we do not know in 
general if the matrices B and C correspond to bounded linear operators and if these operators 
have bounded inverses. See [2] and [3] and references given there for some generalities about 
existence of L[/-factorizations of bounded linear operators as a product of a lower triangular 
and an upper triangular matrix, both corresponding to bounded linear operators. See [3^ f° r 
an example of a unitary operator on ^ 2 (Z>o) without Li7-factorization. 

4 LU factorizations: examples 
4.1 Little g-Jacobi 

From H2.2|) and 1)2, 8|) we see that (|3.13|) will have meaning with the following substitutions: 

y x :=q x , Zl := q~* + abq l ~ x (x,i€Z> ), ^\ := p a ^(x), (4.1) 

PnyUv) 

where 

y v = lim y x = 0, p n {y v ) = (-a)"™ q~H^) ^k. (4.2) 



Rewrite ()3.13|) as 

PniVx) 
Pn{Vu) 



. . min(n,a;) 

E My*))- 1 (nfcta - %)) ^ (nS%) (n£(i - 2/7V)) (4-3) 



fc=0 
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and compare with ()2.2UI) . (|2.21l) . We obtain that 

^ = ife 1 (i-v7 1 w-). ^iu = ^)r 1 (nS(^-'»)) **(n£>w 

Hence, by <223), (SH and (jOl) . 

(b;q)m(.q;q) m ' 



(4.4) 



We can alternatively compute 5 m from (|3.17jl . which takes after substitution of Q4.1)) and 
1)4. 2 jl the form 



(ab,q;q) r 



1 ab, (qab) 2 , — (qab) 2 , a, 0, g' 
(g _1 a&)2, — (g _1 a6)^, 6, a&g r 



;g,g-V" 



(4.5) 



The 6<^4 can be evaluated as a confluent limit case of the summation formula for a terminating 
very well-poised 6^5 series given in (2.4.2)]. The resulting explicit formula for 5 m by 1)4.5)1 
coincides with 1)4.4)) . 

Note that ()3.13)1 with substitution of ()4.1)) and l)4.2|) does not immediately allow to take 
limits for q J, 0. For this we have to renormalize PniVx) by division by Pn{Vv) and we have to 
transfer some factors only depending on fc from Ck,x to B n% k5k / Pn{yv) ■ 

Next we consider formula 1)3.21)) (the upper times lower factorization) with substitution of 
dUD) and (|01> - Then and w n in (fXTTT) become 



a,b;q 



With these substitutions, ()3.21)) can be written in the form 



;(fc-x)(fc-a;-l) 



(-a) 



k~x i ab q n+k ;q)oo(b;q) k f k.-l\ /j.„-i 



/c=max(7i,x) 



(a;g)oo (b;?)^ (g;g) fc 



(^rtte^rt ( 4 - 6 ) 



\{n—x)(n—x—l) ( _„\n— x 



(a^ 2n ;g)oo (b; q) n (g n ; 



71. _— In 



X 2<P2 



L , a6g 

(abq n + x ;q) 00 {q x \q~ 1 

(o;g)oo 



(-a) 

'/" '•V . , 

n-x+1 nhr,2n' q i q a 



(a;g)oo (b; q) x 
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x n x+l 



bq x , q 



,n- n+1 



•n+x ' 



if n > x, 
if n < x. 



(4.7) 



Above we started with ()3.21)) . Instead we might have started with 1)3.18)1 and then obtain 
for C and B ( 1)3.5)1 and ()3.1U)) with substitution of 1)4.1)) ) the inverse matrices by [7J (4.2) and 
(4.11)]. 

Since ()3.21)) was only derived in a formal way, we have not yet proved now 1)4.6)1 and 1)4.7)1 in 
a rigorous way. However, the expressions 1)4.7)1 can be alternatively obtained from 1)2.6)1 by first 
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applying 9 ( (III. 8) or Exercise 1.15 (ii)] and next inverting the series (see Exercise 1.4(H)]): 



{a;q)n V 



■2"i I , _ i i ^ ; ^ 2 "(aft) 1 



g|(n-a:)(n-a:-l) (_ )n-a; 

(<? x ;<? _1 )„ fbq x ,q~ n 



(a;q) n " ri y^-lgl-n-* 



(a;q) n (b;q) x \ « n x+1 



k (a;q)n 
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-,x— n+1 



;q,aq 



(4.8) 



(4.9) 



if n < x. 



Then the expressions (j4.9|) yield expressions (|4.7|) by means of [HI (1.5.4)]. 
4.2 g-Hahn 

See standard formulas for q-Hahn polynomials in 13, §3.6]. They are given by 

]~ n , abq n+1 ,x 



Q n (x;a,b,N;q) := 3 (f> 2 



aq,q 



-N 



(iVeZ >0 , ne{0,l,...,iV}). (4.10) 



It will turn out that we have to renormalize the parameters a and b in order to be able to 
take the limit of these polynomials for q [ 0. For the same reason, we have to consider these 
polynomials with argument q x ~ N , rather than the usual argument q~ x . We define: 



x-N -1 -U at. _\ i I 1 i aoa ii 



-N 



;q,q 



Q% b ' N;q (x) := Q n (q x - I ";q~ L a,q~ 1 b,N;q) = 3 2 

\ «)«/"' 

(N € Z >0 , n, x e {0, 1, . . . , N}; < a < 1, b < 1 or a, b > q 1 ~ N ) 
By [131 (3.6.2)] the functions (|4.11|) satisfy the orthogonality relation 

J2 Q^" 1 ^) Q a n b ' N ' q (x) W a /' N ' q = - " 
x=0 



(4.11) 



(m,n = 0,l,...,iV), 



(4.12) 



where 



IV 



a,b,N;q 



a,b,N;q 



(a;q)N (q N ;q % (b;g) x qX 

(ab;q) N (aq N - 1 ;q- 1 ) x (q;q) x 
(q N ;q- 1 ) n l-abq 2 "- 1 (a,ab;q) n 



(abq N ;q) n 1 - abq^ 1 a n {q,b;q) n ' 



(4.13) 
(4.14) 



Note that the weights w x ' b,q and the dual weights uj% ,q are positive under the constraints for a 
and b given in ()4.11j) . 
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If we apply the transformation formula [HI (3.2.2)] to the terminating 3^2 in (|4.11|) then we 
obtain: 

g -In(n-l) { _ a) -n Q a,b,N;^ x) 
\P'i Q)n 



(q n ,abq n 1 ,q 

= 3< H b,q- N 


-X 

;q, 


qX-N+i\ 
a J 


(4.15) 


( q x , q n , a~ 1 b~ 1 q 


-71+1 




(4.16) 


= C v ~ 1;<rl (^-z) 






(4.17) 


= Q n (q x ;qb-\qa-\N 






(4.18) 


= R x (q n + a- 1 b- l q~ n+ 




-\qa-\N | 1 


7- 1 ), (4.19) 



which is both a polynomial of degree n in q x for n = 0,1, ... ,N and a polynomial of degree 
x in q~ n + a&g™" 1 for x = 0, 1, . . . , N. In KW\ we have used the notation for dual q l -Hahn 
polynomials (see |131 §3.7] for dual g-Hahn polynomials; exchange there in |131 (3.7.1)] n and x, 
replace q by q~ l , and next replace 7 by qb -1 and 8 by (?a _1 ). 
Formula (|4.15j) can be rewritten as: 

n a,b,N;q f ^\ _ \^ } ±(n-ft)(n-fc-l) r_„\n-fc {bq k ;q) n -k (abq n ~ 1 ;q) k (g";^" 1 )^ (q X ]q~ 1 )k 

Qn [ j " ^ 9 1 J (a;g)n(g;g) fc (^KT 1 )* ' 

(4.20) 

In notation (|4,11|) and Q2.2JI the limit formula ^3 (4.6.1)] reads 

lim Qf^(z)=<^(*). (4.21) 

A<^oo 

Expressions (|4. llj) . (|4.15f) . ()4.2U|) for g-Hahn polynomials and expressions (|4.13l) and ()4.14j) for 
their weights and dual weights are very similar to expressions ()2.2|) . ()2.6|) . (|2.9() . (|2.4|) . (|2.5() . 
respectively, for little g-Jacobi polynomials, and the g-Hahn expressions immediately turn into 
their little g-Jacobi counterparts under the limit transition (|4.2H) . 

Remark 4.1. Formula ()4.11|) together with (|4,15|) is the g-Hahn case of a more general identity 
observed by Terwilliger j^Tj related to so-called parameter arrays (see Theorem 4.1 (ii), Lemma 
4.2 and Example 5.4 in [17]). The limit for N — > 00 of this g-Hahn case is the little g-Jacobi 
equality of (|2.2I) and (|2.6|) . This observation may have relevance for the open problem raised in 
|2Tl Problem 11.1]. 

From (|4.11j) and (|4.15|) we see that (|3.13|) will have meaning with the following substitutions: 

y x :=q x , zr^q-^+abq*- 1 (x, i £ {0, 1, . . . , N}, ^jN: := Q% b ' N ' q (x), (4.22) 
where 

Pn(yv) = (-a)- n q--^ n - 1 >^. (4.23) 

(0; q)n 
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With these substitutions, formula ()3.13)) coincides with formula (|4.2U|) if we put 

(qa^r 



(4.24) 



(&; ?)m (?!?)m (9 JV ;9~ 1 )m 

We can alternatively compute <5 m from (|3.17jl . This yields 8 m as the right-hand side of Q4.5JI 
with an additional factor l/(q N ; q~ 1 ) m , from which we obtain (|4.24|) . 

Formula 1)3.21)) can also be specified in the g-Hahn case. Make substitutions as above and, 
furthermore, put 



w x := u n := {{p n {y N )fu% b > N «. 



Then (jHSU) yields 



N 

Q% b > N n(x) = Yj ? 5<*- ai K*-»-i) (-a)* - * 

fc=max(n,x) 



(a;q) N (b;q) x {q;q) k (q N ;q l ) n {q N ;q %' ^ 

Formula l|4.25jl can also be directly reduced to (|4.11|) by first reversing the direction of summation 
in (|4.25l) : substitute k = N — I. Then we obtain 



;q,q 



1 „\ / r,n—N „x—N _ — lu— 1„— N— n+1 

a,b,N;q (x) _ h(N-x)(N-x-l) ( )N~x Igg ? gJiV-x , / g ,9 ' ^ 

(4.26) 

Formula Q4.26JI follows from (|4.11l) by the transformation formula 9, (3.2.2)]. Note that (|4.15f) 
was obtained from 1)4.11)1 by a different application of this transformation formula. Also note 
that the summation reversion changed the upper times lower formula Q4.25JI into the lower times 



upper formula (|4.26j) . 

The s4>2 in Q4.26JI can be written both as a Hahn polynomial and a dual q-Hahn polynomial: 
we rewrite 1)4.26)1 for n, x = 0, 1, . . . , N as 

q -$(N-*XN-*-l) (- a f- N (gigl^zg Q a > b > N «(x) 

{bq x ;q) N -x 

= Q N ^ n (q x - N ; b- l q- N , a~ l q~ N , N; q) (4.27) 

= R N ^ x {q n - N + a- 1 b- 1 q- N - n+1 ;b- 1 q- N , a~ x q- N \ q). (4.28) 

Thus the left-hand side of the above identities is both a polynomial of degree N - x in q n ~ N + 
a~ 1 b~ 1 q~ N ~ n+1 and a polynomial of degree TV — n in q x ~ N . 

Remark 4.2. Analogous to our observation for 1)4.15)1 (see Remark 14. 1)1 . formula 1)4.26)) can also 
be obtained as a consequence of Theorem 4.1 (ii), Lemma 4.2 and Example 5.5 in [57]. Also 
observe that finite orthogonal polynomial systems whose duals are also orthogonal polynomial 
systems, the so-called Leonard pairs, were extensively studied by Terwilliger, see for instance 
)25j . |26j . |27| . |28j . Associated with a Leonard pair is a split decomposition, which gives rise to 
a parameter array. Formula (10) in [27J, which depends on the parameters from that array, is 
essentially the same as our formula (|3.13[) . 
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4.3 O-Hahn functions 

We obtain as an immediate corollary of (|4.2Uj) : 
Theorem 4.3. The limit functions (O-Hahn functions) 



Q% b ' N <°(x) := lim Q% b ' N *(x) (n, x = 0, 1, . . . , N) 



exist. They are equal to the little 0-Jacobi functions pn b, °(x) (see (|2.11l) - (|2.13|) ) restricted to 
x = 0,l,...,N: 



(4.29) 



Q a /' N '°(x) = p a ^\x) (x = 0, 1, . . . , N). (n,x = 0,l,...,N) 
From (J4~T3|) and (}4~l"i|) we get limits 

1-a 



(4.30) 



w 



a,b,N;0 ._ Um w a,b,N;q 
qlO 



l-ab 

(l-a)(l-&) 
l-ab 



a,b,N;0 ._ lim u a,b,N; q 



q\0 



1-b 
l-ab 1 

1 

1 -a 



N 



o(l - b) 
(1 - a)(l - ab) 
a n (l - b) 



if x = 0, 
a 31 if < x < N, 

ifx = N, 

if n = 0, 
if n = 1, 



(4.31) 



(4.32) 



if 2 < n < N. 



Note that 



^o.fc.iViO > o (0 < x < iV) and w^' ff;0 > (0 < n < N) if < a < 1 and b < 1. 
They are almost the same as the weights ()2.18j) and duals weights (j'2. 19f) for little 0-Jacobi 
polynomials. The orthogonality relation (|4.12|) remains valid for q = by taking limits for q [ 0. 



Remark 4.4. Theorem 14.31 was first stated by Haran [111 (7.3.37)], where the limit functions 
(|4.30|) are given in |111 (4.4.10)]. The limit result there follows (although this is not explicitly 
stated) from the expression [111 (7.3.20), (7.3.21)] for g-Hahn polynomials, which reads in our 
notation as: 



Q% b ' N «(x) 



{q N - x ~ n+l ;q)n (b;q) n 



(q N - n+1 ;q)n (a-V~ n ;9)n 3<?>2 V b, qN-x-n+l 



, a q , q 



min(n,x) 

^ q k{n-k)(n 

k=0 



\n-k („N-x. „-l\ , (u„k. 



\q,q 



(4.33) 



+2a;-3fc-l) (~ a )" (g 



(q N ;q 1 ) n (a;q) n _ k (q;q) k 



-{q \q )k(q ;q )*• (4.34) 



Formula (|4.33|) follows from (|4.15f) by the transformation formula 9, (3.2.5)]. Theorem 14.31 can 
be obtained from (|4.34|) by letting q J, 0. 
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Haran's [111 Ch. 7] notation is connected with ours by: 



^ '— = Qf^(x), C W W = ((^l^oo)- 1 ■ (4.35) 



5 Interpretation as spherical functions over j>adic spaces 

We mention for completeness the overall picture concerning our main object of study, the little 
g-Jacobi polynomials. Let F be a local field. F can be Archimedean (R, C) or non- Archimedean, 
that is, either a finite extension of the field Q p of p-adic numbers or the Laurent series over a 
finite field (see |231 Chapter 4] for details). For F non- Archimedean, let O stand for the ring of 
integers. Let be the maximal compact subgroup of GL(d, F): the orthogonal group Od, the 
unitary group Ud or GL(d, O). The natural representation of K-g arising from its action on the 
projective space is given by 

'■ ► U{L 2 ^)), [ P ¥ (g)f](x) = /0T 1 *). 

This representation admits a multiplicity free decomposition into irreducible representations: 

neZ> 

The label F on the various objects here emphasizes the dependence on the field. However, the 
point here is that the decomposition does not depend on the field. Moreover, the irreducibles 
occurring in the decomposition for fixed n correspond to each other when we go through the 
various fields (cf. jllj . |Hj). This correspondence is realized by the observation that, for all F, 
the little g-Jacobi polynomial of degree n has limits which are spherical functions in tl^ for all 
F's. The orthogonality measure of these limit functions is the projection of the Haar measure 
from iff to the space P^ _1 Xk v Fp _1 ; on which the spherical functions live. 

It also turns out (cf. [Oj) that this scheme could be generalized to representations arising 
from the action of these groups on Grassmannians. 



5.1 Interpretation of little O-Jacobi functions 

Let F be a p-adic field, O the ring of integers, p the maximal ideal in 0, p r the cardinality of 
the residue field O/p (p a prime number), and pO = p e (e the ramification index, see again 
|2HI Chapter 4] for details). We look at the representation of GL(d,0), the maximal compact 
subgroup of GL(d,¥), defined by 

p . G L(d,0) — B (L\W d -\¥))) , [p(g)f](x) = /(g^x), 

arising from the action of GL(d,0) on P^^F) = P d-1 (C). Let P m stand for the intersection 
of a standard maximal parabolic subgroup of type (m,d — m) in GL(d,¥) with GL(d,0). In 
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particular, P d — GL(d,0) / P±. When we look at P m -invariants in the representation then 

we have 

L 2 (p d - 1 (0)) Pm = L 2 {P m \GL(d,0)/P 1 ). 

The group GL(d,0) acts on F d ~ l {0) and hence on its quotients F d ~ 1 (0/p k ). Denote the 
stabilizer of (1 : : • • • : 0) G F d - 1 (0/p k ) in GL(d,0) by p[ k) . The space F d ~ 1 (0/p k ) is the 
boundary of a ball of radius k in the rooted tree with root valency p p rZ\ (the- cardinality of the 
projective space over the residue field) and remaining vertices of degree p rd + 1. The orbits of 
P m on F d ~ 1 (0/p k ) consist of k + 1 points and the orbits of the limit space are parameterized 
by Z> U {oo}: 

/ 1 - Q\ 

P m \GL(n, 0)/P x = \\ P m l ° ; ° ";' ° Pi ~ Z> U {oo}. (5.1) 

fc =0 \p k ■■■ 1/ 

The projection fi p of the Haar measure to the orbit space is given (see 111] ) by 

MW) = <' b ''° (k G Z> , a = p-^ d - m \ b = p~ rm ), ^ p ({oo}) = 0, (5.2) 

where the weights w^' , explicitly given by ()2.18j) . are limits for q J, of the weights for the 
orthogonality of the little q-Jacobi polynomials. Thus, one obtains an interpretation of these 
weights on a p-adic space if a, b are as in ()5.2|) . Moreover, for these values of the parameters, the 
p-adic spherical functions (i.e., fixed vectors in the representation p under P m ), are the functions 
Pn' b, ° given by (|2.11j) - ([2.1.3]) . Hence, they are limits of little g-Jacobi polynomials, as was shown 
in [TT]- 

5.2 Interpretation of little O-Laguerre functions 

For the interpretation of the little g-Laguerre polynomials at q = 0, we look at the action of 
GL(m,0) on O m . The orbits of this action are characterized by the minimal valuation of the 
entries of a vector. Hence, the orbits are: O m \ pO m , pO m \ p 2 O m , • • • , {0}. The measure of 

r)~ rrn O'O 

the orbits is precisely the orthogonality measure for the O-Laguerre functions pn ' ' : 

p{pPO m \ pP +1 O m ) = (1 - p-™)p-*3 m (j G Z> ). 
Moreover, by taking GL(m, C)-invariants in the decomposition to irreducibles of the represen- 

Of \ <Q — rm q.q 

tation L (O m ), we find the O-Laguerre functions pn ' ' as fixed vectors. In the special case 
m = r = 1 (group of p-adic units acting on the ring of p-adic integers) this interpretation was 
already obtained by Dunkl and Ramirez jSJ. 

The two pictures (Jacobi and Laguerre) are related in the following manner. One can restrict 
the action of GL(d, O) on P^O), to an action of its subgroup GL(m, O) embedded in the top 
left corner on the subspace {{x : : • • • : : 1)| x G O m } C P d_1 (C). This action is clearly 
the same as the action GL(m,0) on O m . In terms of the parameters a and b, this restriction 
amounts to setting 6 = 0, thus ignoring the irrelevant part of the space. 
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5.3 Product formula — p-adic 

In this subsection we derive the p-adic product formula. We assume that a = p- r ( d - m ) an d 
b = p~ rm . Then the nonnegativity conditions ()2.26[) for the product formula (|2.23l) are valid. 
Let 



"(i) = J2vpU) = { 1 



1 



1 - ab 



if i = 0, 



a 1 if i > 0. 



(5.3) 



Our first step is to look at the spherical functions with a different normalization, which makes 
them idempotents in the convolution algebra L 1 (P rn \GL(d,0)/Pi). We also rewrite them in 
terms of the measure, rather than in terms of a and b. 



a, 6:0/ \ 
e ^ X ) = „ a,6 ;0| | 2 PO ( x ) 



\\P0 



a,b;0 1 \ 
/ \ Pti v'^v a 6*0/ \ 
e n\ x ) = ■ a,6;0|| 2 Pn ' ' ^ 



1, 



(5.4) 



IIP' 




1 



v{n — 1) 
1 1 



u(n) v(n — 1) 
Let {cj = lii,...^}} anci {<?« = }- For i > we have: 



if0<x<n— 1, 
if x = n — 1, 

if x > n — 1 



(n > 1). (5.5) 



Q - Cj + i 
1 

Ci - 



!/(*) 



1/ u 



1) 



Ci-l 



(c_ 



0) 



(5.6) 
(5.7) 



The multiplication in the algebra, which is defined by declaring that the ej's are idempotents, 
is given by: 



d*Cj = v(m&x{i,j})c min{id} , 

f //(max{i,;'})3min{tj} if i ^ j, 

- v(i + l))gi + n(i) J2j>i 9j if « = J- 



(5.8) 
(5.9) 

(5.10) 



At this point we restrict to the p-adic case. In particular we use the fact that the ej's and gi's are 
dual bases in the sense that the former is an idempotent basis for the convolution product, and 
the latter is an idempotent basis for the pointwise product. The spherical transform intertwines 
these products and bases. It follows that the multiplication table for the g^'s (after normalizing) 
with respect to the convolution product is the pointwise multiplication for the idempotents, 
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giving the desired product formula. If we normalize the g^s to be orthonormal, by setting 



Which agrees with the dxyz in (|2.24j) . 

References 

[1] R. Askey and M. E. H. Ismail, A generalization of ultraspherical polynomials, in Studies in 
Pure Mathematics, P. Erdos (ed.), Birkhauser, 1983, pp. 55-78. 

[2] K. T. Andrews, Ph. W. Smith and J. D. Ward, LU -factorization of operators on l\, Proc. 
Amer. Math. Soc. 98 (1968), 247-252. 

[3] K. T. Andrews and J. D. Ward, LU -factorization of order bounded operators on Banach 
sequence spaces, J. Approx. Theory 48 (1986), 169-180. 

[4] R. Askey and J. Wilson, Some basic hypergeometric orthogonal polynomials that generalize 
Jacobi polynomials, Mem. Amer. Math. Soc. (1985), No. 319. 

[5] N. M. Atakishiyev and A. U. Klimyk, On q-orthogonal polynomials, dual to little and big 
q- Jacobi polynomials, J. Math. Anal. Appl. 294 (2004), 246-257; this is a shortened version 
of the preprint arXiv:math.CA/0307250 v3. 

[6] U. Bader and U. Onn, Geometric representations of GL{n, R), cellular Hecke algebras and 
the embedding problem, J. Pure Appl. Algebra (2006), in press; arXiv:math.RT/0411461. 

[7] L. Carlitz, Some inverse relations, Duke Math. J. 40 (1973), 893-901. 

[8] C. F. Dunkl and D. E. Ramirez, A family of countable compact P* -hyper groups, Trans. 
Amer. Math. Soc. 202 (1975), 339-356. 

[9] G. Gasper and M. Rahman, Basic hypergeometric series, Encyclopedia of Mathematics and 
its Applications, Vol. 35, Cambridge University Press, 1990; Second revised edition, Vol. 
96, 2004. 

[10] W. Groenevelt, Bilinear summation formulas from quantum algebra representations, 
Ramanujan J. 8 (2004), 383-416; |arXiv : math . C )A/020 12721 

[11] M. J. S. Haran, The mysteries of the real prime, London Mathematical Society Monographs, 
New Series, Vol. 25, Oxford University Press, 2001. 

[12] H. T. Koelink, Askey-Wilson polynomials and the quantum SU(2) group: survey and ap- 
plications, Acta Appl. Math. 44 (1996), 295-352. 



9i = -tfj)9u w e get: 




if i = j. 



(5.11) 



20 



[13] R. Koekoek and R. F. Swarttouw, The Askey-scheme of hypergeometric orthogonal polyno- 
mials and its q-analogue, Report 98-17, Faculty of Technical Mathematics and Informatics, 
Delft University of Technology, 1998; http://aw.twi.tudelft.nl/" koekoek/askey/. 

[14] T. H. Koornwinder, J acobi functions as limit cases of q-ultraspherical polynomials, J. Math. 
Anal. Appl. 148 (1990), 44-54. 

[15] T. H. Koornwinder and U. Onn, Lower-upper triangular decompositions, q = limits, and 
p-adic interpretations of some q-hypergeometric orthogonal polynomials, 
arX i v : math . RT/ 0405309 v2, 2004. 

[16] C. Krattenthaler, A new matrix inverse, Proc. Amer. Math. Soc. 124 (1996), 47-59. 

[17] I. G. Macdonald, Spherical functions on a group of p-adic type, Ramanujan Institute, 
Madras, 1971. 

[18] I. G. Macdonald, Symmetric functions and Hall polynomials, Oxford University Press, 
Second Edition, 1994. 

[19] I. G. Macdonald, Orthogonal polynomials associated with root systems, 

Seminaire Lotharingien Comb. 45 (2000), article B45a; arXiv: math. QA/001 1046 

[20] A. Okounkov, (Shifted) Macdonald polynomials: q-integral representation and combinatorial 
formula, Compositio Math. 112 (1998), 147-182; [arXiv : q- alg/9605013| 

[21] U. Onn, From p-adic to real Grassmannians via the quantum, Adv. Math. 204 (2006), 
152-175; larXiv : math . RT/04 05138, 

[22] G. Polya and G. Szego, Aufgaben und Lehrsatze aus der Analysis, zweiter Band, Dritter 
Auflage, Springer 1964; English translation, Problems and theorems in analysis. II, Springer, 
1998. 

[23] D. Ramakrishnan and R. J. Valenza, Fourier analysis on number fields, Graduate Texts in 
Mathematics 186, Springer- Verlag, 1999. 

[24] H. Rosengren, A new quantum algebraic interpretation of the Askey-Wilson polynomials, 
in: q-Series from a contemporary perspective, M. E. H. Ismail and D. W. Stanton (eds.), 
Contemporary Math. 254 (2000), 371-394. 

[25] P. Terwilliger, Two linear transformations each tridiagonal with respect to an eigenbasis of 
the other, J. Linear Algebra Appl. 330 (2001), 149-203. 

[26] P. Terwilliger, Leonard pairs from 24 points of view, Rocky Mountain J. Math. 32 (2002), 
827-888; larXiv : math~R A/0406577l 

[27] P. Terwilliger, Two linear transformations each tridiagonal with respect to an eigenbasis 
of the other; comments on the parameter array, Des. Codes Cryptogr. 34 (2005), 307-332; 
arXiv:math. RA/0306291, 



21 



[28] P. Terwilliger, Leonard pairs and the q-Racah polynomials, Linear Algebra Appl. 387 (2004), 
235-276; |arXiv : matn~C) A/030630 1 

[29] N. Ja. Vilenkin and A. U. Klimyk, Representation of Lie groups and special functions, Vols. 
1, 2, 3, Kluwer, 1991, 1993, 1992. 

[30] M. Voit, A product formula for orthogonal polynomials associated with infinite distance- 
transitive graphs, J. Approx. Theory 120 (2003), 337-354. 

[31] R. A. Zuidwijk, Complementary triangular forms for infinite matrices, in Operator theory 
and boundary eigenvalue problems (Vienna, 1993), Oper. Theory Adv. Appl. 80, Birkhauser, 
1995, pp. 289-299. 



Tom H. Koornwinder 

Korteweg-de Vries Institute, University of Amsterdam, 

Plantage Muidergracht 24, 1018 TV Amsterdam, The Netherlands 

thkOscience . uva . nl 
Uri Onn 

Einstein Institute of Mathematics, Edmond Safra Campus, Givat Ram, 
The Hebrew University of Jerusalem, Jerusalem 91904, Israel 

urionOmath . hu j i . ac . il 



22 



